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We discuss the stability of homonuclear and heteronuclear mixtures of 3 He and 4 He atoms in the 
metastable 2 3 Si state (He*) and predict positions and widths of Feshbach resonances by using 
the Asymptotic Bound-state Model (ABM). All calculations are performed without fit parameters, 
using ab-initio calculations of molecular potentials. One promising very broad Feshbach resonance 
(AB = 72.9iig - 3 mT) is found that allows for tuning of the inter-isotope scattering length. 

PACS numbers: 



I. INTRODUCTION 

The helium atom is one of the most simple atoms. Its 
electronic structure with only two electrons allows ab- 
initio calculations of level energies with extreme preci- 
sion, testing basic theory of atomic structure. Also for 
the interaction between two helium atoms ab-initio quan- 
tum chemistry calculations allow highly accurate molec- 
ular potentials in some cases. In this paper we focus 
on helium atoms in the metastable ls2s 3 Si state (He*) 
for which molecular potentials have recently been calcu- 
lated [THl] . We investigate the possibilities to modify col- 
lision properties of mixtures of He* atoms in magnetic or 
optical dipole traps. Bose-Einstein Condensation (BEC) 
of 4 He* atoms has been realized by several groups [5]- 
[9]. This metastable isotope of helium has no nuclear 
spin, and is magnetically trappable in the fully stretched 
\s,m 3 ) = |1,+1) state, where s and m s are the electronic 
spin and magnetic quantum numbers respectively. The 
fermionic isotope 3 He*, which, due to its nuclear spin 
z=l/2, shows hyperfine structure, has been cooled in the 
f=3/2 hyperfine manifold (with f = s + i) to degeneracy 
by sympathetic cooling with 4 He* [TO] . 

In the ultracold regime the scattering length a accu- 
rately describes the interaction between the atoms. Mag- 
netically tunable Feshbach resonances [TT] can be utilized 
to, in principle, tune the scattering length at will. Nu- 
merous experiments have been proposed based on this 
tuning possibility. For 3 He*- 4 He* mixtures it would be 
possible to observe phase separation when the scattering 
lengths are large and positive [12 . Using a position sen- 
sitive micro-channel plate detector to detect single atoms 
(with close to 100% efficiency) may reveal boson-induced 
p-wave pairing of the 3 He* fermions [13]. This detection 
technique has also allowed an atom-optics detection of 
the Hanbury Brown and Twiss effect for both bosons [14] 
and fermions [15 , and the possibility to tune interactions 
may also open up new research possibilities within this 
field PI]. 

All experiments with ultracold He* atoms so far have 
been performed in fully stretched low-field seeking states 
(|1,+1) for 4 He*, \f,m f ) = |3/2,+3/2> for 3 He*). The 
atom-atom interactions for 4 He* in the fully stretched 



spin states are known with astonishingly high preci- 
sion [4j [17] . From the binding energy of the v = 14 ro- 
vibrational state in the quintet potential, which is only 
slightly below the dissociation threshold, a large positive 
scattering length a=7.567(24) nm was deduced [il fTUlfTT] . 
in very good agreement with experimental findings [18]. 
For collisions between atoms which are not in the fully 
stretched spin states, Penning Ionization (PI) losses will 
strongly compromise the stability of trapped He* atoms 
in experiments. 

For the above mentioned reasons it has been favor- 
able experimentally to prepare 4 He* atoms in the fully 
stretched spin state. However, the absence of a nuclear 
spin in 4 He* strongly limits the possibilities to tune the 
inter- atomic interactions via Feshbach resonances. Al- 
though it is possible to induce Feshbach resonances via 
the magnetic dipole interaction, a more efficient coupling 
to other states occurs when an internal hyperfine struc- 
ture is present in at least one of the two interacting 
atoms. This occurs when we allow for mixtures between 
4 He* and 3 He*. A recent experiment has already shown 
that 4 He* atoms can be trapped in an optical dipole 
trap [19J, potentially allowing trapping of all magnetic 
substates of both helium isotopes. 

In this paper we discuss the possibility to access Fesh- 
bach resonances in collisions between 4 He* and 3 He* 
atoms. In addition, we also will consider homonuclear 
3 He*- 3 He* and 4 He*- 4 He* collisions. We will work with 
the Asymptotic Bound-state Model (ABM) [20] which 
will be reviewed in SectjH] To determine the binding 
energies of the molecular potentials for the three possi- 
ble isotope combinations, we use the accumulated phase 
method [21 . These binding energies will serve as in- 
put parameters for the ABM. In order to discuss exper- 
imentally relevant systems, we will evaluate the various 



possible loss mechanisms in Sect fIV] after which, possible 
Feshbach resonances will be explored in SectfV} 



II. THE ASYMPTOTIC BOUND-STATE MODEL 

To predict the magnetic field position and width of 
Feshbach resonances, we use the Asymptotic Bound-state 
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Model (ABM) [20] [22j [23]. We recapitulate the basic 
principles of the ABM and show how this model can be 
applied to metastable helium. For an elaborate discus- 
sion on the ABM we refer the reader to Ref. [20] . In our 
approach, we ignore the effect of PI on the elastic inter- 
actions. A more detailed discussion of PI follows later 
in Sec. |IV| To simplify the notation in discussing various 
helium isotope combinations, we use that f = s in the ab- 
sence of nuclear spin for 4 He*, thereby allowing similar 
notational treatment of the spin structure of fermionic 
and bosonic helium atoms. 

The ABM enables us to determine the energy of the 
coupled molecular states, the eigenstates of the total 
Hamiltonian H, without solving the actual coupled radial 
Schrodinger equation. For the collision of two metastable 
helium atoms in an external magnetic field the total 
Hamiltonian is 

H= l^+Hint + V + Vdd, (1) 

where the first term represents the relative kinetic energy 
with /i the reduced mass, and Hint the two-body internal 
energy. The internal energy of 3 He* and 4 He* as a func- 
tion external magnetic field is shown in Fig. [I] Note the 
inverted hyperfine structure of 3 He*, which means that 
the / = 3/2 manifold is below the / = 1/2 manifold. 
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Figure 1: Single-atom internal energy diagrams for 3 He* and 
4 He*. The capital (lower case) letters are used to label the 
energies of 4 He* ( 3 He*). The magnetic substates are also 
labeled by their /, mf quantum numbers, where at non-zero 
magnetic field m/ is still a good quantum number, whereas / 
is not. 

The central (Coulomb) interaction V depends on the 
magnitude of the total electron spin S = Si + S2 and the 
distance between the nuclei r = |r|. This interaction can 
be decomposed as 

2 

V(r) = ^ 25+ V(r)P 5 , (2) 
s=o 

where 25+1 V(r) is the adiabatic molecular potential for 
the molecular state 25+1 £+ n , and Ps projects onto the 



S'-subspace. For 4 He*, short range molecular potentials 
have been calculated ab-initio by Miiller et al pQ. In 
this paper, we refer to the molecular potentials as singlet 
S = ( X S+ ), triplet S = 1 ( 3 S+ ) and quintet S = 2 
( 5 S+ ) potentials. This nomenclature is based on the 
spin configuration of all four electrons of the two He* 
atoms. Note that often a single He* atom is also referred 
to as being in a spin-triplet state. Here, however, we 
do not use this designation to avoid confusion with the 
molecular singlet /triplet /quintet potentials. 

The direct dipole-dipole interactions of the electronic 
spins [47 Vdd can be written as a scalar product of two 
irreducible spherical tensors of rank 2: 

3a 2 +2 

V dd (r) = -— MHr ® r} 2 _ g { Sl s 2 } 2 „ (3) 

T q=-2 

where {r (g) r}2™ = yff ^2m(0 is a spherical har- 
monic [24], and a is the fine structure constant. Being 
anisotropic, the interaction allows for redistribution of 
angular momentum between spin S and orbital £ angu- 
lar momentum. Therefore, although much weaker, this 
interaction can couple many more states as compared to 
isotropic interactions. The orbital angular momentum 
coupling can only occur for states when t — t' = 0, ±2, 
with the exception of £ = £' — which is forbidden. The 
spin angular momentum coupling by Vdd obeys similar 
selection rules; S — S f = 0, ±2, where S = S f = is not 
allowed. 

From ultracold collisions it is well known that the po- 
sition of the least bound state is crucial for the determi- 
nation of the interaction properties. In the ABM we fol- 
low the same philosophy, and since the vibrational level 
splitting of the least bound levels for a light atom such 
as metastable helium is much larger (at least h x 21 GHz 
for the 4 He*- 4 He* system, where h is Planck's constant) 
than the hyperfine energy, we only use the properties 
of the least bound state. We expand the coupled so- 
lutions into a product basis of spin states and molecu- 
lar states of the uncoupled 25+1 V(r) potentials. These 
molecular states are the eigenstates of the relative Hamil- 
tonian, defined as H re \ = p 2 /2/i + V, and are denoted 
as l^^) with corresponding binding energies ef £ . We 
have defined A G {33,34,44} to label the different iso- 
tope combinations which we consider. Having specified 
the Hamiltonian and the appropriate basis, the energies 
of the coupled bound states follow from a simple ma- 
trix diagonalization as a function of the strength of the 
magnetic field. 

To determine the characteristic properties of Feshbach 
resonances, i.e. their magnetic field widths and positions, 
we need to examine the behavior of the coupled bound 
states near the threshold of an open channel. A channel 
is called open when the total energy exceeds the chan- 
nel threshold energy, which is defined as the sum of the 
single atom internal energies. The intersection of the 
coupled bound state with this threshold gives the posi- 
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tion of a Feshbach resonance, which is accurate to the 
order of the width of this resonance. However, in reality 
the coupled bound state will acquire an increasing open 
channel component as it approaches the open channel 
threshold. Therefore, near threshold the binding energy 
curve will bend quadratically towards the threshold curve 
as a function of magnetic field [11]. From this behavior 
of the bound state it is possible to extract the width, and 
an improved estimate of the position of the resonance. 

Since we distinguish between open and closed chan- 
nels, it is useful to partition the space of states describ- 
ing spatial and spin degrees of freedom into an open and 
a closed channel subspace [25j [26] . The Hamiltonian for 
the system is written as 



T~L — Hpp + Hqq + Hpq + 7~Lqp, 



(4) 



where U PP = PUP.Upq = PUQ, etc. Here P and Q 
are projection operators onto the open and closed chan- 
nel subspace respectively. The bound states of Hqq will 
be responsible for Feshbach resonances when coupled to 
the open channel subspace via HpqCHqp)- These bare 
Q-space bound states, denoted by |0q), then become 
dressed by this interaction. Usually this dressing oc- 
curs by coupling to scattering states in the open channel 
P-space whereas for the ABM we use the bound states 
\<pp) of the open channel subspace. The magnetic field 
at which the energy of the dressed bound state becomes 
degenerate with the energy of the threshold is where a 
Feshbach resonance will occur. 

The width of the resonance will depend on the cou- 
pling strength between the open and the resonant closed 
channels (Hpq), and the binding energy of the open- 
and closed-channel bound states (ep,6Q respectively). 
We define the width AB as the difference in magnetic 
field where a = and a = oo. For this purpose we 
introduce the ^-matrix, which can be written for elas- 
tic scattering as S = e 2lS ^ where S(k) is the scatter- 
ing phase, with Hk = y/2jlE, and the energy E is the 
collision energy defined with respect to the open chan- 
nel threshold energy. The scattering length is defined as 
a = — lim/c-^o tan5(fc)/fe. 

The energy of the dressed bound state, which is also 
referred to as the Feshbach molecular state, corresponds 
to a pole of the total ^-matrix. This scattering matrix 
S = SpSq (of the effective problem in P-space) is a prod- 
uct of the Sp matrix, and a resonant part Sq which in- 
volves coupling of a Q-space bound state to P-space [27] . 
The Sp scattering matrix, which describes the scattering 
process in P-space in the absence of coupling to the Q- 
space, is determined by considering only the dominant 
bound state in P-space [2D] . Here we neglect the other 
nearby resonance poles in P-space [28] which, is a valid 
approximation if the background scattering length (ab g ) 
is larger than the typical range of the interaction poten- 
tial [20]. By determining the total S- matrix, we are able 
to deduce the (magnetic field) position and width of the 
Feshbach resonance. 



The Sq scattering matrix is usually determined by a 
single Q-space bound state. However, it is possible that 
multiple Q-space bound states have to be taken into ac- 
count to properly describe a Feshbach resonance, as will 
be the case for the wide 3 He*- 4 He* resonance discussed 
later in Sect. El 

The multiple Q-state expression for Sq is given by [26] 



S Q = 1- 2m 
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(n) 



i ~ E — En 



(5) 



where E n are the complex eigenvalues of Hqq + Wqq. 
The operator Wqq = Hqp E +^ Hpp Hpq describes a 
temporary transition from Q-space to P-space, prop- 
agation in P-space, and re-emission into Q-space [48] . 
These eigenvalues can be determined by solving the sec- 
ular equation in the basis of Q-space bound states \(/)q). 
By demanding unitarity of the ^-matrix the coupling el- 
ements 7^ n ) can be expressed in terms of the complex 
energies E n [26]. Hereby we have completely specified 
the total ^-matrix. 

The free parameters in our model which determine the 
magnetic field resonance position Bo and the field width 
AB, are the binding energies Cg £ and the overlap be- 
tween various molecular states (^^l^/^/) [20]. Next 
we discuss how we can obtain these quantities by utiliz- 
ing known molecular potentials. 



III. MOLECULAR STATES 

The essential parameters for the ABM are the binding 
energies ef £ of the molecular potentials. If these values 
are known, Feshbach resonances can be predicted with 
an accuracy determined by the accuracy at which the 
ef £ parameters are known. Vice-versa it is possible to 
obtain the binding energies by fitting the calculated res- 
onance positions to experimentally observed resonances 
as has been shown in Ref. [22]. For metastable Helium 
no Feshbach resonances have been observed so far. To 
predict resonance positions we thus require knowledge of 
the binding energies of the S = 0, 1,2 potentials. These 
potentials are known from literature, however there is a 
significant difference between the S = and S = 1 po- 
tentials on one side and the S — 2 potential on the other 
side. The former ones are usually described as complex 
potentials to incorporate the effect of PI and have been 
described to a certain degree of accuracy by Miiller, et 
al. [I]. The latter potential (S = 2) does not exhibit 
PI and can therefore be accurately described by a purely 
real potential, which has been measured and calculated 
with very high accuracy [4j [18] . 

In order to properly account for these uncertainties 
we make use of the accumulated phase method [21] [29] . 
Moreover, since we study different isotope combinations, 
the uncertainties of the potentials are reliably treated 
via mass scaling within this method. The s-wave accu- 
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mulated phase can be written as 

A 4>s(E) = A <pf(E) + A A 4> s , (6) 

where A (j)^\E) is the uncorrected accumulated phase, 
resulting directly from the radial wavefunction in the in- 
ner region of a molecular potential 2S + 1 V ? for a given 
energy E, at a radial range tq. We allow for a phase cor- 
rection A A (j)s, which accounts for the mismatch between 
calculated and experimentally measured quantities. It 
is determined from an asymptotic boundary condition, 
by demanding a particular energy for the highest bound 
state in the potential. The accumulated phase for a dif- 
ferent isotope combination A' is found by a mass-scaling 
of the phase correction: 

A '4>s(4') = A '<t> { s\4') + KA A <Ps, (7) 

where 1Z = \J \±a> / \±a- We have verified the accumulated 
phase calculations by comparing with calculations per- 
formed on the full (complex) potentials. In the following 
we will discuss the different potentials in more detail. 

S=2 potential 

A highly accurate ab-initio S = 2 potential was de- 
termined by Przybytek and Jeziorski [4j[l7]. For 4 He* - 
4 He* they predicted a binding energy of the least bound 
state = 14) equal to -89.6(8) MHz. Using two- 
photon photo-association spectroscopy, Moal et al. [18] 
measured the binding energy of this least bound state to 
be er, 4 /h = —91.35(6) MHz. We use this measurement as 
input parameter for 4 He* - 4 He* ABM calculations, and 
we use Eq. ^ to find the binding energies for the other 
isotope combinations. 

We construct the full S = 2 potential by fitting the 
short-range potential by Przybytek [4j to the long-range 
dispersive potential V^i sp (r) = — C^/r 6 — Cg/r 8 — Cio/r 10 
using the accurately known dispersive coefficients [30] . 
We forge these potentials around r = 20 ao by vertically 
shifting the short-range potential to match the long-range 
potential and by applying a smoothing function. By the 
demand of a bound state at e^/h = —91.35(6) MHz, 
we apply the phase correction at ro = 18 ao, which is 
determined with Eq. The results for the differ- 

ent isotope combinations are e^/h = —4.84 MHz and 
e^/h = —413.83 MHz. These energies compare well to 
those obtained by Przybytek and Jeziorski [31] based on 
the full scattering potentials. 

S=0,1 potentials 

The short range parts of the S = 0,1 potentials, as ob- 
tained by Miiller et al. [1 , are known accurately enough 
to calculate the position of the least bound levels in these 
potentials with relatively large errorbars. The S = 0, 1 
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Figure 2: Variation in the binding energies of the least bound 
state of the S = 0, 1 potentials. The uncertainty in e^ 4 
leads to a range of mass-scaled binding energies, indicated 
by thick black lines. Within the present accuracy, the triplet 
3 He* - 3 He* potential can become deep enough to capture a 
new bound state. This weakly bound state cannot be dis- 
tinguished on this scale. The inset demonstrates how we 
find the binding energies for this case. The accumulated 
phase of Eq. ^ (which is indicated by a grey shaded area) 
is matched with the phase of the wavefunction (solid line), 
starting from an asymptotic bound state boundary condi- 
tion. These accumulated phases are plotted as a function 
of kf = y/\el 3 /h\ (MHz 1/2 ). The variation in the accu- 
mulated phase, caused by the uncertainties of the poten- 
tial, corresponds to a range of binding energies. Therefore, 
kf 3 can take any of the values < kf 3 < 1 MHz 1/2 or 
175 < kf < 215 MHz 1/2 

potentials also include an imaginary part to incorporate 
for the PI losses. As the ABM works with only real bind- 
ing energies we assume real S = 0, 1 potentials. This in- 
troduces an additional error to which we will come back 
later. The use of the accumulated phase method allows us 
to incorporate these different uncertainties in the A A <fis 
parameter. 

We use results obtained by Leo et al. [32J to set up- 
per and lower bounds on the scattering lengths for these 
potentials 

af = 34(10)a , a 44 = 32(9)a , (8) 

where the error bars indicate the maximum inaccuracy 
of both the real and complex short range potentials as 
found by Leo et al. [32] . The scattering lengths enable us 
to determine binding energies for the bosonic metastable 
helium system, as well as for the other isotope combina- 
tions. The results of these calculations are summarized 
in Fig. [2] 

As mentioned before scattering in the S = 0, 1 poten- 
tials experience strong PI losses which are generally de- 
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scribed by including an imaginary term in the potential. 
Up till now we have neglected this imaginary terms, how- 
ever, to calculate the binding energies we use the scatter- 
ing lengths (Eq. ^) which are based on the full optical 
potentials. Therefore, this procedure accounts partly for 
the inaccuracy in the binding energies induced by using 
a real scattering potential. 

Additionally, we have verified that using only real po- 
tentials induces a relatively small error by comparing the 
value of ef 4 obtained with the accumulated phase method 
with the value obtained using the complete S = 1 poten- 
tial. For the imaginary potential we use the autoioniza- 
tion width as given in Ref. pQ. The total S = 0, 1 po- 
tentials were constructed by vertically shifting the short 
range potentials pQ to (smoothly) match 5 V(r) — V exc h(r) 
for r > 11.5 ao- The exchange terms, which depend on 
quantum number 5, were determined in Ref. [33 J. We 
find that including the autoionization width shifts the 
real part of e| 4 by ~ h x 1.5 GHz. This is a significant 
shift, but still smaller than the uncertaincies in the bind- 
ing energies caused by the errors as given in Eq. 
Therefore, the limiting factor in the present calculations 
is the inaccurately known short range potentials rather 
than the approximation of using only real potentials. 

For the S = 0, 1 potentials of 3 He*- 3 He* we find that 
the scattering lengths are negative; i.e. virtual bound 
states will dominate the low energy scattering properties 
of these potentials [28] . Within the variation of the S = 
1 potential, it is even possible that the virtual bound 
state turns into a (weakly) bound state, and therefore 
becomes the highest bound state of the potential which 
will make the scattering length (large) positive. This is 
displayed in the inset of Fig.|2j The low energy scattering 
properties of the S = 1 potential will depend strongly on 
the precise position of such a weakly bound state, which 
will inhibit accurate predictions for Feshbach resonances 
without more accurate knowledge of the S = 1 potential. 

IV. INELASTIC DECAY PROCESSES 

To determine which two-body hyperfine states are best 
suited to explore Feshbach resonances, we discuss various 
loss mechanisms that may occur. A particular combi- 
nation of states can be relatively stable, however, the 
occurrence of Feshbach resonances may limit the stabil- 
ity, since the underlying resonant bound state can suffer 
strongly from inelastic effects. We discuss the possibility 
of resonantly enhanced losses in the next section. 

In contrast to ground- state atoms, metastable helium 
atoms can undergo highly exothermic ionizing collisions, 
since the internal energy of two He* atoms exceeds the 
He* ionization potential by 15 eV. These reactions, de- 
scribed by 

He * +He ^{S4 + H e e -, + e " ® 
will be referred to as Penning Ionization (PI), which in- 



cludes the process commonly referred to as associative 
ionization. As this reaction is electrostatic, the total 
electron spin S is conserved (Wigner spin-conservation 
rule). For fully stretched states, where 5 = 2, the re- 
action would violate spin conservation and is therefore 
forbidden in first order, whereas when S = or S = 1 
the reaction can proceed. The probability for PI for these 
latter potentials is ~0.975 [1 , hence a severe loss process 
for collisions involving scattering through the singlet or 
triplet potential. 

Another important loss mechanism is spin exchange 
relaxation which is induced by the central part V of the 
interaction. For these isotropic interactions the projec- 
tion mp = rrif 1 + m/ 2 of total spin angular momentum 
F = fi + f2 on the magnetic field axis is conserved dur- 
ing the collision. By preparing atoms in the energetically 
lowest two-body hyperfine state, for a particular value of 
m]?, only unfavorable endothermic collisions can occur, 
effectively suppressing spin exchange relaxation losses. 
In Fig. [3] we have labeled these two body states by their 
one-body constituents. 
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Figure 3: The two-body hyperfine diagrams are shown for 
3 He*- 4 He* and 3 He*- 3 He*. For each possible value of mf, 
the energetically lowest two-body hyperfine state (i.e. stable 
against spin exchange relaxation) is labeled (black), whereas 
the other states (unstable against spin exchange relaxation) 
are not labeled (gray). For homonuclear fermionic helium, 
symmetry prevents that hyperfine states composed of identi- 
cal substates (dot-dashed) can be populated with t — even. 

The spin-dipole interaction Vdd between the spins 
of both electrons induces losses as well. For these 
anisotropic interactions the projection of total angular 
momentum T = F + i is conserved, coupling only two- 
body states of the same mjr = + vng. Here, two 
mechanisms may cause trap loss: spin relaxation (a re i) 
and relaxation-induced ionization (a r i). The first pro- 
cess only relaxes the spin projection m^, whereas for the 
second process the S — 2 state is coupled to S = 0, 
which decays via the 'normal' PI mechanism. At low 
temperatures and low B fields the latter process dom- 
inates. These mechanisms were found to be the most 
prominent cause of losses in a spin-polarized gas of 4 He*. 
The loss rate was calculated to be four orders of mag- 



6 



nitude smaller [34], as compared to an unpolarized gas, 
where direct PI is dominant for trap loss. 

Generically three-body loss rates depend on the value 
of the two-body scattering length. For a homonuclear 
gas, where 3 (or 2) identical bosons or 2 identical fermions 
participate in a three-body collision, the three body loss 
rate will vary as |a| p , where p > 3 [35 . Therefore large 
values of \a\ are predicted to (strongly) enhance three- 
body loss rates, which are affected by spin state and 
statistics of the participating atoms in the process. 

The dominant loss mechanisms for trapped metastable 
helium atoms are Penning ionization and spin exchange 
relaxation. Therefore we will only consider two-body 
hyperfine states which are stable against these decay 
processes. Since symmetry will impose additional con- 
straints, we will elaborate on the loss mechanisms for the 
homonuclear and heteronuclear case separately, in order 
to find the right experimental conditions (a sufficiently 
long lifetime) to search for Feshbach resonances. Future 
experiments based on predictions of the ABM discussed 
in Sect.[ll]will depend heavily on these background losses. 

Homonuclear losses 



the gas against PI losses at fiK temperatures. For 3 He*, 
studies of losses in a MOT have also been performed and 
also here theory and experiment agree on the loss rate 
in the dark: i^3 npol) =2xl0- 10 cm 3 /s [37 . In the dark 
a mixture of all four magnetic substates of the trapped 
f=3/2 hyperfine manifold shows a loss rate that is even 
larger than for 4 He*. The loss rate here is described very 
well assuming only ionizing collisions between atoms in 
different magnetic substates. 

The nuclear spin of 3 He* gives rise to magnetic field 
dependent PI loss. Therefore, we have plotted the S = 2 
fraction of the two-body hyperfine states (stable against 
spin exchange relaxation) as a function of the magnetic 
field in Fig. |4j This decomposition into molecular states 
containing quantum number S of the two-body hyperfine 
state allows us to estimate the Penning ionization rate 
which will be large for states containing a small S = 2 
fraction [37] . For s- wave collisions of two fermions, the 
a + b state {mp = — 2) will scatter via dominantly the 
5 = 2 potential at higher magnetic fields whereas for all 
the other states unfavorable PI processes will dominate. 
If we consider £ =odd collisions as well, the two-body 
states a + a, d + d are stable against PI losses. 



For homonuclear collisions (between identical bosons 
or fermions), the symmetrization requirement 

# + / + £ = even, (10) 

for their two-body state can have severe consequences 
on the stability of the gas. If we for example consider 
collisions of 4 He* atoms, Penning ionization losses via 
the 5 = 1 potential can only occur via £ =odd colli- 
sions. These losses are therefore strongly suppressed at 
T ~ l^K, where the experiments around degeneracy take 
place [36] . 

For bosonic helium atoms all two-body hyperfine states 
except the degenerate B + B,A + C states are stable 
against spin exchange relaxation. Additionally, PI will 
also induce losses for these two states making them un- 
suitable in our pursue of Feshbach resonances. The fully 
stretched states A + A, C + C however, are stable against 
PI losses. If we restrict the temperatures such that 
we can consider only s- wave collisions, the A + B, B + C 
states will also only scatter via the 5 = 2 potential. The 
stability of this mixture will however be limited by the 
losses between atoms in the B state. The stability can 
be improved by making B the minority spin species. In 
a magneto optical trap PI losses have been studied in 
the presence of MOT light and without. Good agree- 
ment between theory and experiment was obtained for 
loss rates of unpolarized atoms in the dark; the two-body 
loss rate turned out to be very large: K^ npo1 ^ =lx 10 -10 
cm 3 /s [37]. 

For fermionic helium atoms, ultracold collisions be- 
tween atoms in the same substates can only occur via 
odd partial waves (Pauli principle), effectively stabilizing 



In the experiments on degenerate gases of either 3 He* 
or 4 He*, three-body losses have been shown to be less 
important than two-body losses. In the case of the fully 
stretched C + C state a three-body loss rate constant 
L=2xl0 -27 cm 6 /s 2 was calculated [34], which experi- 
mentally was confirmed in BEC lifetime studies [7] for 
4 He*. 



3 He*- 4 He* 3 He* - 3 He* 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

B(T) B(T) 



Figure 4: The quintet fraction for two-body hyperfine states 
as a function of magnetic field is shown for 3 He*- 4 He* and 
3 He*- 3 He*. We only show states which are stable against 
spin exchange relaxation, i.e. the energetically lowest two- 
body state for a particular tuf value. The a + a and d + d 
states (gray) can only have collisions with £ =odd due to 
symmetry. 



7 



Heteronuclear losses 

For collisions between 3 He* and 4 He* atoms we do not 
have symmetrization requirements as for the homonu- 
clear case. Only for atoms in the fully stretched states 
PI is suppressed [37] [38]. In a MOT containing both 
isotopes a heteronuclear ionization rate coefficient was 
deduced from loss measurements that also agrees with 
theory: i^4 npol) =3x KT 10 cm 3 /s [38 . This shows that 
losses in an unpolarized He* gas at mK temperatures, for 
all isotopes, are well understood and may be extrapolated 
to jaK temperatures. 

Considering heteronuclear collisions, we have to take 
into account that only 3 He* and 4 He* substates should be 
selected which are intrinsically stable. The previous dis- 
cussion on the homonuclear situation then dictates that 
4 He* atoms in the B state have to be excluded. The 
decomposition into the quintet fraction of the possible 
entrance states as a function of magnetic field is shown 
in Fig. [4] We find that in addition to the fully stretched 
a + A, d + C states, the b + A state will contain a large 
quintet fraction. 

As the quintet scattering length is extremely large, we 
expected that the stability of a boson-fermion mixture 
of spin-polarized 3 He* and 4 He* atoms is severely com- 
promised by three-body recombination processes. Ini- 
tial observations indeed confirm this [10]. Recent exper- 
iments on a boson-fermion mixture of 87 Rb- 40 K Simoni, 
et al. [39 however, have shown that such a system can 
be made more stable by having an excess of the fermions 
rather than the bosons. This is explained by the fact 
that for a three-body loss process in that case two iden- 
tical fermions have to come close to each other which is 
suppressed by the Pauli principle. 

Effect on Feshbach resonances 

In view of inelastic loss processes, the selected sta- 
ble two-body hyperfine states for all three helium iso- 
tope combinations are: 4 He*- 4 He*: {A + A, C + C}, 
3 He*- 3 He*: {a + a, a + b, d + d}, and 3 He*- 4 He*: 
{a + A, b + A, d + C}. Before we discuss the possibility 
of finding Feshbach resonances for these two-body states, 
we will consider the impact that inelastic loss processes 
might have on these resonances. 

The presence of inelastic loss processes such as spin re- 
laxation, relaxation-induced ionization, and Penning ion- 
ization will affect the Feshbach resonances. These inelas- 
tic events may occur not only in the open V but also 
in the closed Q channel subspace. The scattering length 
describing the two-body interactions, including these ef- 
fects, will now be complex valued [40], and the divergence 
of the real part of the scattering length at resonance will 
be suppressed. The strength of the resonance will depend 
on the relative magnitudes of the coupling elements be- 
tween the resonant state to the elastic and inelastic chan- 
nels [41]. 



Although Feshbach resonances are usually associated 
with various enhanced (two and three body) loss pro- 
cesses, they can also have a stabilizing effect [42] [43] . 
Since inelastic losses can be induced in both the V- and 
Q-space by PI, it is possible that PI losses can be sup- 
pressed as the admixture of V- and Q-space can change 
in the vicinity of a Feshbach resonance. 

For PI losses, two metastable helium atoms (in close 
proximity of each other, for S ^ 2) couple to energeti- 
cally lower ionic states which yields an inelastic process. 
To describe the effect of PI losses, one usually takes op- 
tical potentials and avoids the use of the ionic channels. 
This transforms the closed quantum system to an open 
one, i.e. the effective Hamiltonian describing the collision 
between two helium atoms has become non-Hermitian. 
This is different as compared to the inelastic scattering 
caused by the magnetic dipole-dipole interaction, which 
is described by increasing the number of channels in the 
open channel subspace but keeping the effective Hamil- 
tonian Hermitian. Although the description of both pro- 
cesses is different, the effect of these processes will (in 
principle) be the same. 

The PI process will also influence the molecular states 
of the S = 0, 1 potentials. The generic effect of this 
imaginary potential to describe PI, is that it will cause 
the bound states of the real potential to become unsta- 
ble [44], i.e. they acquire a finite lifetime. The binding 
energies of such unstable bound states are complex val- 
ued, as one would expect for a non-Hermitian Hamilto- 
nian. The imaginary part of the complex energy of such 
an unstable bound state is not used as a parameter in 
the ABM. 



V. FESHBACH RESONANCES 

Considering the various inelastic decay mechanisms 
due to PI and spin exchange relaxation, we have nar- 
rowed down the number of interesting open (also referred 
to as entrance) channels dramatically. In our search 
for Feshbach resonances our focus will be on resonances 
caused by coupled bound states which are mainly in an 
S = 2 state. The reason for this lies in the inaccuracy 
of the S = 0, 1 potentials which will lead to a significant 
spread of possible singlet and triplet binding energies [49J , 
as has been discussed in Sect. IIIIl We limit the search for 
Feshbach resonances to magnetic fields up to 1 T. At the 
end of this section, the found Feshbach resonances are 
summarized in Table [I] 

Feshbach resonances induced by the magnetic dipole- 
dipole interaction Vdd will a ls° be considered, although 
these are expected to yield much weaker Feshbach reso- 
nances as compared to the ones induced by central inter- 
actions V. Since the dipole-dipole interaction is weak we 
consider only first order processes and thus only include 
partial waves £ < 2. The basis set for the ABM calcu- 
lations will consist of the (bound) eigenstates and 
energies e^ £ of the 25+1 V(r) + £{£ + l)h 2 /(2/iir 2 ) poten- 



8 



tials for £ < 2. These molecular states are determined in 



a similar fashion as presented in Sect III for s- wave bound 
states. The eigenstates IV^) allow us to determine the 
matrix elements of Vdd as 

(*£i\v dd \H,i>) = -3a 2 (^i^i^,^) E (-i) q 

1 q=~2 

'-{lm g \Y 2 ^ q {r)\i f m t ){(i\{^ ® s 2 } 2q W),(ll) 



where |^^, £f ) = \ipg' i )\£m£)\a), and \a) is the two-body 
spin state [20]. The r-dependent factor (ipf i\r~ 3 \ijj°s f £/ ) 
is determined by a numerical integration. By diagonaliz- 
ing the Hamiltonian H we find the energies of the coupled 
bound states, as described in Ref. [20]. The magnetic 
dipole-dipole interaction is thus treated as a perturba- 
tion in first order in the ABM, similar to the treatment 
of magnetic dipole-dipole interactions for metastable he- 
lium atoms by Beams et al [45] . 

We selected our entrance channels by minimizing pos- 
sible inelastic losses due to PI and spin exchange relax- 
ation, however, Vdd niay couple to other states which 
are significantly less stable than these entrance chan- 
nels. In view of the Feshbach formalism discussed in 
Sect. [II] these inelastic processes occur when the dimen- 
sion of the open channel subspace dim(V) > 1. The 
inelastic losses induced by Vdd will alter the Feshbach 
resonance characteristics, e.g. the real part of the scat- 
tering length will not diverge on resonance as it would 
if there was only elastic scattering [41 . Throughout this 
section, just as with bare closed channel bound states, 
we only use the dominant (i.e. energetically closest to 
the threshold) open channel bound state for the determi- 
nation of Bq,AB. We neglect the effect of the inelastic 
losses due to Vdd and PI on the Feshbach resonances. 



Homonuclear gas 

For 4 He*- 4 He* collisions, the absence of nuclear spin 
prohibits (coupled) bound states within an m^-manifold 
to cross the open channel threshold and induce Feshbach 
resonances. Since Vdd can couple states of different m^, 
Feshbach resonances can be induced by the spin-spin in- 
teraction. For the two-body state A + A we find two 
Feshbach resonances. It is important to note that the 
entrance channel A + A is purely quintet and there is a 
bound state e^/h = —91 MHz parallel to the entrance 
channel. 

A d-wave singlet bound state couples the entrance 
channel via Vdd and causes a narrow Feshbach resonance 
at Bq = 9.9 mT whose (field) width equals AB = 0.2 /aT 
which is denoted as I in Fig. [5] If we make the singlet 
potential deeper, within the inaccuracy of the potential, 
we may find this resonance for fields up to Bq = 46 mT 
with the same field width. By making the potential more 
shallow, in comparision with the nominal potential, the 



resonance will be found at lower magnetic fields and can 
even dissappear. Around B = 133.1 mT a p-wave triplet 
bound state will cross the threshold of A + A (denoted 
by II). This bound state can however not couple to the 
open channels and will therefore not cause a Feshbach 
resonance. For higher magnetic fields, a Feshbach reso- 
nance is found at Bq = 546.0 mT and a field width of 
AB = 1.1 jaT. This resonance is caused by a d-wave quin- 
tet bound state which is shown at III in Fig. [5] Since the 
S = 2 potential is the most accurate potential we have, 
we expect that the prediction of the III resonance will be 
the most accurate, as opposed to I where an S = bound 
state causes the resonance. For the C + C channel, Fesh- 
bach resonances are absent. We note however that within 
the uncertainty of the singlet potential it is possible for 
a S = 0,£ = 2 bound state to be relatively close to the 
S = 2,£ = bound state at low magnetic fields. Since 
these two bound states can be coupled by the Vdd inter- 
action, the S = 2,£ = bound state will become less 
stable. This can be of interest to recent studies [45j [46] 
on the lifetime of the S = 2,£ = bound state of the 
C + C channel. This qualitative argument needs further 
investigation. 

For two fermionic helium atoms colliding, the en- 
trance channels a + a, d + d can only be populated for 
£ =odd. Therefore, for ultracold scattering experi- 
ments, the only relevant s- wave channel is a + b. For 
this (anti-symmetrized) state the scattering is domi- 
nated by the quintet potential at high magnetic fields, 
see Fig. [4] Without threshold effects, we expect that 
around B « 970 mT an s-wave quintet bound state 
will cross the scattering threshold. Unlike the Feshbach 
resonances discussed so-far, the energy of this coupled 
bound state is higher than the threshold energy for low 
magnetic fields. Including threshold effects we find that 
this bound state will not cause a Feshbach resonance. 
There is an s-wave triplet bound state which does cre- 
ate a wide (AB ^15.2 mT) Feshbach resonance around 
Bq ~ 1426.5 mT. The magnetic dipole-dipole interaction 
induces five Feshbach resonances which are all caused by 
d-wave bound states in either the singlet or triplet con- 
figuration, hence these resonance positions will not be 
stated here. 

Since the S = 1 potential is almost resonant within 
the uncertainty variations, it may be able to capture a 
new bound state, see the inset of Fig. [2] This makes it 
very challenging to reliably predict Feshbach resonances 
for 3 He*- 3 He* that involve the 5=1 potential. 



Heteronuclear gas 

For the heteronuclear gas the only selected channel 
which allows for Feshbach resonances induced by the cen- 
tral interaction is b + A. The a + A and d + C mixtures 
are both fully-stretched states, and therefore can only 
have Feshbach resonances induced by Vdd- If we apply 
the ABM without taking into account threshold effects, 
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Figure 5: (Color online) Energies of coupled bound states 
(black curve) are plotted versus the external magnetic field 
for 4 He*- 4 He*, neglecting threshold effects. The energy of the 
A + A, A + B states correspond to different channel threshold 
energies (gray curve). At the points I, II, III a coupled bound 
state intersects with the A + A threshold. The threshold of 
A + A is degenerate with a d-wave A + A state. For the points 
I, III Feshbach resonances are induced by d-wave bound states, 
with 5 = 0,2 spin quantum numbers respectively. For I there 
will be a large spread in possible resonance field position Bo 
because of the inaccuracy of the singlet potential. At II a p- 
wave bound state crosses the threshold but this state cannot 
induce a Feshbach resonance as it is not coupled the entrance 
channel. 



we find two resonances: one at low field Bq = 1.3 mT and 
one Bo = 347.8 mT. However, since the least bound quin- 
tet state (e^ 4 ) is almost resonant, threshold effects will 
dominate and broad Feshbach resonances are expected. 

To include threshold effects we apply the theory as 
described in Sect. [TTJ We determine the uncoupled V, Q 
space bound states. Remarkably we do not find that 
either of the bare closed channel bound states (with en- 
ergies 6q x and €q 2 ) cross the scattering threshold at low 
magnetic field (see Fig. [6] ), as one would usually ex- 
pect. Counter intuitively it is the energy of the bare 
open channel bound state ep that becomes degenerate 
with the scattering threshold which, if we could physi- 
cally uncouple V- and Q-space (Hpq — >• 0), would re- 
sult in a potential resonance. This will have severe con- 
sequences for the observed resonance structure. Where 
usually a Q-space bound state pushes the dressed bound 
state through threshold, here it is the P-space bound 
state. 

The single resonance approximation will now fail 
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Figure 6: (Color online) In the upper panel the scattering 
length is shown for the b + A state as a function of magnetic 
field. The scattering length is scaled to the van der Waals 
range r v dw = 34 ao- At a magnetic field of B = 121.4 mT 
the Feshbach resonance will occur, which has a field width of 
AB = 72.9 mT. At a magnetic field of Bq « 214.7 mT, just 
outside the plot, the 7^-space bound state will become reso- 
nant. In the lower panel the energies of the uncoupled V, Q- 
space bound states are shown (dashed lines blue ep and red 
curves €q 1 ,€q 2 respectively) with respect to the open chan- 
nel threshold energy. The energy of the dressed bound state, 
which causes the Feshbach resonance, is found by solving the 
pole equation of the ^-matrix (black solid curve). 



since the |0q x ) and \4>q 2 ) states are energetically al- 
most equidistant to the threshold. To describe the 
dressed bound state of the coupled system we need to 
study the peculiar interplay of three bare bound states: 
|0p)j \4 > Q2)- T ne interplay of these bound states is 

illustrated nicely by the pole equation of the total scat- 
tering matrix S = SpSq: 



(K P + ik)(E-E 1 )(E-E 2 ) = 0, 



(12) 



where we have used Eq. ([5| for Sq. The closed channel 
bound states |0q 1 ), \(j)Q 2 ) not only interact with the \<pp) 
state but also with each other via P-space. The energy of 
the dressed bound state which results from this interplay 
is shown as a solid line in Fig. [6| Here only the physical 
solution of Eq. ( 12 ) which causes the resonance is shown. 



From these threshold effects, we find a Feshbach res- 
onance position Bq 



121.4^45.9 mT and width AB 
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72.9l}g"3 mT. The variation in the position and width 
of the resonance stated here is due to the uncertainty 
in the triplet bound state, since a significant fraction of 
the |0q 1 ) state is in a triplet state. The strong coupling 
between open and closed channel states (large %pq) in 
combination with the P-space bound state being almost 
resonant, results in this very broad Feshbach resonance. 

The open channel bound state crosses the threshold at 
a mag netic field of = 214.7 mT. The description of 
the binding energies and the scattering length for B > 
Bq will become inaccurate. 
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Figure 7: (Color online) Coupled bound states diagram for 
the 3 He*- 4 He* mixture in the a + A, b + A states. Only the 
positions of Feshbach resonances induced by S = 2 bound 
states are labeled and the precise values are given in the text. 
For the spin stretched entrance channel a + A, we find two 
resonances induced by (d-wave) quintet bound states at (1,2). 
For the b + A entrance channel quintet (d-wave) bound states 
cause Feshbach resonances at I, II, III. For IV a p-wave bound 
state crosses the threshold. This bound state can only be 
coupled to other states of ^=odd, and will thus not couple 
to the threshold s-wave state. All other crossings correspond 
to coupled bound states of S < 2. Note the large number of 
possible resonances at low magnetic fields. 

The magnetic dipole-dipole interaction for the b + A 
state can induce multiple Feshbach resonances. Since 
the a + A channel is energetically open when we include 
Vdd we expect the Feshbach resonances to be modified by 
these open channels. Therefore we will only predict the 
positions (thus without threshold effects) of the three res- 
onances Bq = {361.8,572.9,587.9} mT induced by quin- 
tet d-wave bound states, I, II, III respectively in Fig. [7| 
At IV a p-wave quintet bound state will cross the thres- 
hold. This bound state cannot couple to s-wave scatter- 



Table I: A summary of the predicted Feshbach resonances for 
He*. The accuracy of the predicted resonance field position 
(-Bo) is mainly determined by the dominant S value of the 
responsible coupled bound state. All error bars are based on 
the deviations in the used parameters, not on the errors made 
by applying the ABM. 1 ) For A = 44 the resonance at B = 
9.9 mT can, within its error bars, shift up to Bo — 46 mT 
and below Bo = mT, i.e. disappear. 2 ) For A = 33 the 
nominal result is stated, although, due to the resonant nature 
of the S = 1 potential, the resonance structure can change 
dramatically. This will lead to a large deviation (« 300 mT) 
in B . 



A 


mixture 


Bo (mT) 


AB (mT) t S coupling 


44 


A+A 


9.9 1 


0.2xl0~ 3 2 V dd 


44 


A+A 


546.0(1) 


lxlO- 3 2 2 Vdd 


33 


a+b 


1426.5 2 


13.2 1 V 



34 b+A 121.4tg-7 72.9±i|-J 2 V 

34 b+A 361.8(6) - 2 2 Vdd 

34 b+A 572.9(1) - 2 2 V dd 

34 b+A 587.9(1) - 2 2 V dd 

34 a+A 503.0(1) 5.5xl0" 3 2 2 V dd 

34 a+A 994.1(1) O.lxlO -3 2 2 V dd 



ing states in b + A. 

In the fully stretched a + A, d + C states, only Vdd 
induced Feshbach resonances can occur. For a + A 
we find five Feshbach resonances. Two of those are 
caused by quintet (d-wave) bound states which yield 
B = 503.0 mT, AB = 5.5 /iT and B = 994.1 mT, 
AB = 0.1 /iT labeled as 1, 2 respectively in Fig. [7) The 
d + C entrance channel state is expected to suffer more 
from inelastic losses due to Vdd as there are multiple (non- 
degenerate) open channels. We only find one Feshbach 
resonance for this channel which is caused by a triplet 
(i-wave bound state. 



VI. CONCLUDING REMARKS 

We presented here the first study of Feshbach reso- 
nances for metastable helium atoms, by using the Asymp- 
totic Bound-state Model. By analyzing the various in- 
elastic decay processes we have selected a few two-body 
spin states suitable for observing Feshbach resonances. 
Reliably predicting these resonances is in some cases hin- 
dered by fact that the S = 0, 1 potentials are known 
with far less accuracy as compared to the S = 2 molecu- 
lar potential. Therefore, we have limited our discussion 
to coupled bound states which have a dominant quintet 
character and cause Feshbach resonances for magnetic 
fields up to 1 T. To study these resonances we have uti- 
lized and expanded the ABM: magnetic dipole-dipole in- 
teractions as well as overlapping resonances can now be 
described with this model. Although we found several 
Feshbach resonances, we did not find wide resonances for 
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the homonuclear (bosonic and fermionic) gas for the se- 
lected spin states. The heteronuclear system however, 
reveals a very wide resonance AB = 72.9 mT at rela- 
tively low magnetic field Bq = 121.4 mT making it of 
potential interest for further theoretical and experimen- 
tal investigation. 

Measurements of Feshbach resonances will aid enor- 
mously in constructing more accurate S = 0, 1 molecular 
potentials. Vice- versa, with more accurate S = 0, 1 inter- 
action potentials we will be able to give reliable predic- 
tions of the abundant number of resonances induced by 
coupled bound states in a dominantly singlet or triplet 
state. It is our current understanding that these poten- 
tials can induce resonances at lower magnetic fields as 
compared to S = 2 dominated bound states which makes 



them of great potential interest. Accurate interaction 
potentials would pave the way for full, numerically ex- 
act, coupled channels calculations which will yield more 
accurate predictions. Based on qualitative arguments, 
we also point out the interesting possibility of Feshbach 
resonance induced stabilization of PI losses, where the 
resonance can effectively reduce the inelastic loss rate. 
This effect may for instance be used to stabilize S = 0, 1 
(coupled) bound states. 
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